Structure of magnetic order in Pauli limited unconventional superconductors 
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We analyse the magnetic structure of the antiferromagnetic order induced in Pauli limited d-wave 
superconductors by Zeeman coupling to magnetic field. We determine the phase diagram in the H-T 
plane, and find that the magnetic phase, which is stabilized at low temperatures and just below the 
upper critical field, can have two realizations depending primarily on the shape of the underlying 
Fermi surface. The double-Q magnetic ordering may persist over the entire coexistence range. 
Alternatively, there may exist a weak first order transition from a double-Q structure at lower fields 
to a single-Q modulation at higher fields. Together with the calculations of the NMR line-shape 
these results suggest the second scenario as a serious candidate for describing the superconducting 
state of CeCoIns. 
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Introduction. Understanding of the emergence and 
stability of competing orders in correlated systems has 
been a major focus of research. Static long-range or- 
der changes the excitation spectrum of itinerant systems, 
and the resulting energy gain determines which instabil- 
ity dominates. A single order parameter appears in sim- 
ple cases, but electron-electron interactions can also favor 
several distinct phases in many materials. The compe- 
tition between different ordering phenomena make these 
phases sensitive to applied pressure, chemical doping and 
magnetic field. 

Heavy fermion CeCoIns presents one of the most 
prominent and puzzling examples of such complex behav- 
ior. At ambient pressure it is a very clean singlet d-wave 
superconductor with lines of nodes in the gap function. 
Upon doping, superconductivity coexists with and then is 
pre-empted by an antiferromagnetic (AFM) order. The 
isostructural CeRhlns is an AFM metal. At low tem- 
peratures and fields, H, just below the Pauli limited up- 
per critical field H C 2, CeCoIns enters a thermodynamic 
phase that was initially conjectured pQ to be the first real- 
ization of the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) 
state, where the superconducting (SC) order parameter 
oscillates in real space with a length scale proportional to 
the ratio between H and the Fermi velocity. However, ex- 
periments showed that in this superconductivity coexists 
with static long-range incommensurate AFM order [2H7] , 
even though no AFM order is found upon suppression of 
SC by doping, magnetic field, or pressure [8UTT]. Incom- 
mensurability and small ordered moment are consistent 
with itinerant, or spin-density wave (SDW) magnetism, 
which normally competes for the electronic states with 
superconductivity. The inescapable conclusion that in 
CeCoIn 5 superconductivity enables antiferromagnetism 
challenged our views. 

The existing theories fall broadly into two categories. 
The first assumes that the spatial variation of the SC or- 
der under applied field (due to the FFLO modulation) 
yields a higher single-particle density of states (DOS) in 



the regions where superconductivity is suppressed, and 
nucleates a SDW [T2HT5] . However, single NMR line 
shape indicates homogeneous magnetism [5j[16], and, to- 
gether with independence of the SDW wave vector on the 
field [6 , suggests that the origin of the SDW instability is 
not related to this modulation (although the FFLO state 
may still exist [16]). The second category investigates 
how a strong Pauli limiting in a nodal superconductor 
may promote a uniform SDW state [T7lf2T] . We recently 
showed [21] that a two-dimensional d-wave superconduc- 
tor under a Zeeman (paramagnetic) field is generically 
unstable towards formation of an SDW at the wave vec- 
tors ±Q± connecting the opposite nodes of the SC order 
parameter. The instability is due to nearly perfect nest- 
ing of field-induced pockets of Bogoliubov quasiparticles, 
and explains the incommensurate and field-independent 
SDW wave vector, as well as the direction of the local mo- 
ment normal both to the applied field and to the struc- 
tural layers. 

The main challenge to these theories is the structure 
of the SDW state. Very generally, if the instability origi- 
nates from magnetic scattering of quasiparticles between 
the nodal regions, the staggered magnetization, ttiq, of 
a d-wave superconductor has Fourier components for the 
four wave vectors ±Q± connecting the "nested" pairs of 
nodal points, i.e. along each of the two orthogonal di- 
rections, 110 and 110. We term this structure double-Q 
(2Q). The 2Q SDW gaps all four pockets of Bogoliubov 
quasiparticles thus leading to the greatest energy gain. 
Therefore, it is always more stable for weak-coupling 
and hence lower fields [21]. However, the NMR line- 
shape [5j [16] is consistent solely with a single-Q (1Q) 
SDW modulation. We conjectured previously that inter- 
mediate or strong-coupling may stabilize 1Q state, but, 
to our knowledge, no resolution of this discrepancy, which 
is generic to the theories of SDW instability of the nodal 
superconductors, has been offered. 

Below we show that interference between Bogoliubov 
quasiparticles from the neighborhood of different pock- 



2 



ets is the relevant factor that determines the magnetic 
structure. The phase diagram contains a transition line 
between 2Q and 1Q magnetic orders for a range of pa- 
rameters in the strong-coupling limit. This transition has 
weak thermodynamic signatures (such as the specific heat 
anomaly) but manifests itself very clearly in the change 
of the NMR lineshape. Unambiguous observation of such 
a change would strongly favor the current theory for the 
emergence of the SDW order in CeCoIn 5 . 

Model: We consider a mean-field Hamiltonian for a d- 
wave superconductor under Zeeman magnetic field, 
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Here e k — 2t (cos k x + cos k y ) + At' cos k x cos k y — \i is the 
band energy, h = gfisH/2 is the Zeeman splitting , N is 
the number of lattice sites, and A k = A (cos k x — cos k y ) 
is the superconducting order parameter. J denotes mag- 
netic interaction, and we explicitly separated the four 
(equal in magnitude) wave- vectors ±Q V : v = ± in- 
dicates a different pair of "nested" gap nodes, while 
Q v and Q v = —Q v are the two vectors connecting 
a given pair, see inset of Figjlja). Below we select 
Q± — (77r/6, ±77r/6) by fixing the chemical potential 
to \i = £q + /2- The self-consistency equations are 
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FIG. 1: Thermodynamic phase diagrams for (a) t' /t = 0, 
V/t = 3, J/t = 4, and (b) t' /t = 0.2, V/t = 0.8, J/t = 2.4. 
Open and closed circles represent first and second order phase 
transition, respectively. Crosses represent the transition from 
single-Q or double-Q to uniform superconductivity for L = 
432. The inset of (a) shows the nesting vectors Q± and the 
Fermi surfaces at t' /t — (dashed line) and t' /t — 0.2 (solid 
line), with V/t = 0, and J/t = 0. 



where (• • • ) denotes the thermodynamic average for T~L, 
and V is the pairing interaction in the d-wave channel. 

The four order parameter components tuq lead to the 
following real space distribution of magnetic moments 
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The easy axis in CeCoIn 5 is perpendicular to the lay- 
ers [22], so the field in the plane clearly favors the 
staggered moment along the easy direction, in agree- 
ment with experiment [22j [23] . We therefore assume 



such collinear spin ordering, \tuq u \ = \rn-Q |, and choose 
the easy axis as the spin x-direction (crystal c-axis), 
m Qu ~ m h - There is a zero phase mode for each value 
of v due to the incommensurate nature of the ordered 



phase. Without loss of generality we choose the ampli- 
tudes mQ to be real numbers, i.e., we fix both phases at 
zero. Two possible magnetic structures can be stabilized 
under these conditions. The 1Q structure, in which only 
one of the two amplitudes, rriQ u1 is finite, corresponds 
to a SDW that is modulated along the direction paral- 
lel to the corresponding Q v vector. The 2Q structure , 
\ m Q + \ — \ m Q-\, is modulated along both 110 and 110. 

By numerically solving Eqs. Q-([2| on a square lattice 
of linear size L ~ 864 and for different values of band pa- 
rameters and coupling constants, we found two distinct 
types of phase diagrams that are shown in Figs. [2] Both 
of them show the 2Q ordering at low temperatures and 
moderately high fields. However, while this phase per- 
sists throughout the whole magnetic region for some pa- 
rameter values, Fig. [TJa), in other cases the SDW phase 
exhibits a transition between 1Q and 2Q structures, see 
Fig. [TJb). The transition between uniform SC and nor- 
mal (N) phases is the same in both situations: second 



3 



order at low fields/high temperatures, and first order at 
high fields/low temperatures; this behavior is generic for 
two-dimensional paramagnetically limited superconduc- 
tors. The transition from the normal to the magnetically 
ordered SC phase remains of first order, while the onset 
of SDW in the SC phase occurs via a second order tran- 
sition. Finally, a weak first order phase transition line 
separates the 1Q and 2Q phases. Since neutron scatter- 
ing [23] and NMR measurements [5] of CeCoIn 5 clearly 
indicate the presence of a 1Q phase at high fields, this 
latter case is relevant [2T] , and we now discuss the origin 
and the properties of the 1Q/2Q transition in detail. 

Recall that in the weak-coupling (low field) regime, 
the SDW order arises from nesting between opposite 
pockets of Bogoliubov quasiparticles (eigenstates of Ho + 
^bcs) [2] generated by the Zeeman term in a nodal su- 
perconductor [24 . The stability of the 2Q phase arises 
from the decoupling of the mean field equations ([2| for 
the four amplitudes raQ, which requires the same am- 
plitude |t71q| of all the components to fully gap the four 
nodal pockets. Numerical solutions for all sets of param- 
eters always yield the 2Q order at the lowest fields, in 
agreement with this argument. 

The lowest order interaction term of a Ginzburg- 
Landau expansion of the free energy is \rriQ + \ 2 \rriQ_ | 2 . 
Diagrammatically, the coefficient of this term involves 
a product of four propagators of the Bogoliubov quasi- 
particles at momenta separated by ±Q±, with the main 
contribution when all four momenta are near the Fermi 
surface pockets. Since the ordering wave vectors Q± are 
incommensurate, the near-nodal Bogoliubov quasiparti- 
cles do not satisfy this requirement and give a small con- 
tribution to this coefficient. In contrast, once the pockets 
grow to include areas away from the nodes, near (±7r, 0) 
and (0, ±7r) points, the constraint is satisfied. Therefore 
the coefficient of the \rriQ + \ 2 \rriQ_ | 2 is small at low fields, 
but, as the field increases and the pockets grow, the two 
components of the magnetization begin to interfere de- 
structively, leading to an increase in the energy of the 
2Q phase relative to that of the 1Q phase. To test this 
argument we computed the average occupation number 
of Bogoliubov quasiparticles with a given wave- vector k 
over the entire Brillouin zone for the parameters where 
2Q/1Q transition exists. The results, shown in Fig. |2j 
clearly indicate that the 1Q phase becomes stable as soon 
as the overlap between particle-hole clouds from differ- 
ent pockets becomes strong. Therefore, the main physical 
parameter that decides between the two possible scenar- 
ios shown in Figs. [TJa) and (b) is the barrier between 
orthogonal pockets which, in turn, depends on the shape 
of the Fermi surface. Recall that the SDW phase appears 
for values of the magnetic interaction J > 0.85 — 0.9J C , 
where J c is the critical value when the normal state be- 
comes unstable towards magnetic order [21] : this result 
remains valid for any Fermi surface shape. We verified 
that, for a Fermi surface where 2Q-1Q transition exists, 



the 1Q phase persists over wide range of values of the 
pairing interaction V, although, once V is significantly 
reduced, the 1Q region becomes too narrow to be found 
numerically in our calculation. At the same time, for the 
Fermi surface where only 2Q order is found, such as that 
used in the main panel of Fig. [TJa), no changes in the 
interaction can induce the 1Q phase. 
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FIG. 2: Density of Bogolyubov quasi-particles (7^7*.^) for 
t'/t = 0.2, V/t = 0.8, J/t = 2.4, and T/t = 10~ 4 with L = 
864. Shaded panels (c), (d), and (h) correspond to the lowest 
free energy solution of the self-consistent equations ([2| at a 
given magnetic field value. 

Given that the 1Q and 2Q phases can only exist to- 
gether in the phase diagram, we analyse the phase tran- 
sition between the two. Figures [3^a-c) show the so- 
lutions of Eqs.Q-Q at low T for the parameters in 
Fig. [lfb). While both the SC and the SDW order pa- 
rameters show a significant jump at the transition from 
the normal paramagnetic phase, the change in Ao at 
the 1Q-2Q transition is extremely small, and the re- 
duction in rriQ + is compensated by the concomitant in- 
crease in rriQ_ , so that the spatially- averaged magnetiza- 
2 

tion \m r \ = (1/L d )^2 r m 2 = (rriQ + + wiq_)/2, remains 
nearly unchanged. This is also clear from the almost un- 
noticeable kink in the free energy density at the transition 
line. As a result, the latent heat at the transition between 
N and 1Q phases is much larger than the latent heat at 
the 1Q-2Q transition (see Fig. [3jc)). In other words, the 
1Q-2Q transition, while first order, has very weak ther- 
modynamic signatures, and is difficult to detect by bulk 
measurements. On the other hand, probes sensitive to 
the local magnetic structure, such as neutron scattering 
or NMR, should be able to clearly distinguish these two 
different phases. Figure |3|d) shows the predicted NMR 
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FIG. 3: Mean-field solutions for t'/t = 0.2, V/t = 0.8, 
J/t = 2.4, and T/t = 0.003 with L = 432. (a) Order pa- 
rameters (Ao, rriQ + , and rriQ_) and the spatially-averaged 
magnetization \m r \. (b) Free energy densities /mf/£ f° r dif- 
ferent solutions of self-consistent equations (c) Entropy s 
as a function of magnetic field, (d) Expected NMR spectrum. 



line shape in the magnetically ordered Q-SDW phase and 
for different values of the magnetic field, obtained by cal- 
culating the distribution of local magnetic field generated 
by the sum of the four Ce spins that surround each In(l) 
site [5]. There is a qualitative difference between the 
NMR line-shapes for the 1Q and 2Q phases: while the 
1Q phase leads to a double-horn NMR line shape char- 
acteristic of ID modulations, the line shape of the 2Q 
phase has a maximum at the center (logarithmic Van- 
Hove singularity) typical of 2D modulations. The NMR 
line shape as a function of increasing field, see Fig. [3jd), 
is very consistent with recent measurements from two 
different groups [5j [16] . 

Conclusions. In summary, we presented a theory for 
the magnetic structure of the low-temperature and high 
field phase of Pauli limited nodal superconductors such 
as CeCoIns. Our main finding is that such materials 
are expected to broadly fall into two classes depending 
predominantly on the shape of the underlying Fermi sur- 
face: those exhibiting only the phase in which the local 
moments are modulated along the two orthogonal direc- 
tions connecting opposite nodes of the superconducting 
gap, and those exhibiting a transition between such a 
doubly modulated phase and SDW order along a sin- 
gle inter-nodal direction. The origin of the transition 
is in the destructive interference of scattering processes 
of Bogoliubov quasiparticles with two orthogonal wave- 
vectors. Thermodynamic signatures of the 2Q-1Q tran- 
sition are very weak, and therefore local magnetic probes 
are best suited for probing it. The NMR line shape in 
our analysis is in a good agreement with recent data. It 



would be highly desirable to further test the proposed 
transition using neutron scattering. 
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